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;> ■ ABSTRACT 

' LMC X-l and LMC X-3 are the only known persistent stellar-mass black hole candi- 

OO . dates that have almost always shown spectra that are dominated by a soft, thermal 

component. We present here results from 170ksec long Rossi X-ray Timing Explorer 
!■/"") ' (RXTE) observations of these objects, taken in 1996 December, where their spectra 

, can be described by a disc black body plus an additional soft (T ~ 2.8) high-energy 

power-law (detected up to energies of 50keV in LMC X-3). These observations, as 
well as archival Advanced Satellite for Cosmology and Astrophysics (ASCA) observa- 
tions, constrain any narrow Fe line present in the spectra to have an equivalent width 
is 90 eV. Stronger, broad lines (w 150 eV EW, a w 1 keV) are permitted. We also study 
the variability of LMC X-l. Its X-ray power spectral density (PSD) is approximately 
^ ■ oc / _1 between 10~ 3 and 0.3 Hz with a root mean square (rms) variability of w 7%. At 

^ | energies > 5keV the PSD shows evidence of a break at / > 0.2 Hz, possibly indicating 

an outer disc radius of ^ 1000 GM/c 2 in this likely wind-fed system. Furthermore, 
the coherence function 7 2 (/), a measure of the degree of linear correlation, between 
variability in the > 5 keV band and variablity in the lower energy bands is extremely 
low (^ 50%). We discuss the implications of these observations for the mechanisms 
' that might be producing the soft and hard X-rays in these systems. 
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1 INTRODUCTION 



The soft state has bee n observed in steady sources such 



Since t 



1965), 



as LMC X-l and LMC X-3 ( Treves ct al. 198£ - L Ebisawa et al 



198£ ; Treves et al. 199C ; Ebisawa et al. 1993; Bchmidtke et al 



Vio gtndy of o-alartir Mart- hclo ranrHA atca mTTPg) 



has shotvn that these objects display a large variety of states 



199!.), in r ecurring transients such as GX 339—4 (Grebenev 
et al 19951), and in a number of transients such as Nova Mus- 



which are characterized by their distinct spectral shapes and 
temporal behaviours. The most important states which have 
been identified are the "low/hard state", which is character- 
ized by p hard Y-ray gpprtmm with a phi-ytnn indpy T — 1 7 

and lar^S root mean square (rms) variability <~. 30% (Tanaka 
& Lewin 1995; Nowak 1995; and references therein), and the 
"high/soft state", which is spectrally softer (T ~ 2.5) and 
exhibits less variability. The fractional Eddington luminos- 
ity of sources in the soft state tends to be higher than that 
for sources in the hard state (Nowak 1995, and references 
therein). 



cae (Miyamoto et al. 1994). The persistent BHC Cygnus X-l 
has been observed to switch between the hard and the soft 
state (albeit with kT ~ 0.3 keV for the soft state) , with the 
total X-ray luminosity s taying roughly constant (|Cui et al 



Zhang et al. 1997). 



199J ; 

A great deal of observational attention has been focused 
on the more commonly observed hard state, since most of 
the brighter galactic BHCs only occasionally transit to the 
soft state. Only two of the persistent nearby BHCs, LMC X- 
1 and LMC X-3, until recently have a lways been observed 
in the soft state. Wilms et al. (1999d), hereafter Paper II, 
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present evidence that LMC X-3 periodically transits into 
the low/hard state. In this work we present 170ksec long 
Rossi X-ray Timing Explorer (RXTE) observations of both 
LMC X-l and LMC X-3 during high/soft X-ray states. 

LMC X-3 is a highly v ariable BHC with a probable 
9Mq compact object mass (Cowley et al. 1983). Its lumi- 



nosity has been observed to be as high as 4 x 10 38 erg s 1 , 
which i s ~ 30% of its Eddington luminosity, and its soft X- 



ray flux (~ 1-9 keV) is variable by a factor of more than four 
on long time scales (see Paper II). LMC X-3 also has exhib- 
i ted a strong 99- or 198-day periodicity in its soft X-ray flux 
( Cowley et al. 1991 ; |Cowley et al. 1994 ) , which is also appar- 
ent in the ASM monitoring, albeit with a period that varies 
over time (Paper II). Previously we had suggested that this 
periodicity mig ht be associated wi th a warped, precessing 
accretion disc (|Wilms et al. 1999b| ); however, with the re- 
cently observed "state changes" to a low/hard flux it now 
seems likely that a systematic variation of the accretion rate 
is an important part of this long-term variability. 

LMC X-l is also a good candidate for a black hole. 
Using a l arge n umber of ROSAT HRI observations, Cow- 
ley et al. (1995) were able to ident ify th e counterpart with 
"star number 32" of Cowley et al. (1978). This object has a 
mass function of only / = 0.144 M©, but including other 
evidence th e mass of the compac t object appears to be 
M > 4M a (|Hutchings et al. 1987| ), and probably ^ 6M 
(Cowley et al. 1995). The lu minosity of the ob ject is typ- 



ically about 2 x 10 as erg s _1 (Long et al. 1981). Although 



small differences between LMC X-l and LMC X-3 are appar- 
ent, their spectra and short-time temporal behaviours have 
historically been quite similar. LMC X-l, however, does not 
exhibit any obvious periodic behaviour in its long term X- 
ray lightcurve. 

We have monitored LMC X-l and LMC X-3 with RXTE 
in three to four weekly intervals since the end of 1996 in order 
to enable a systematic study of the soft state. The campaign 
started with 170ksec long observations of both sources, and 
in this paper we present results from the spectral and tem- 
poral analysis of these long observations. In addition, we 
consider spectral results of archival Advanced Satellite for 
Cosmology and Astrophysics (ASCA) observations. Prelimi- 
nary results of our analyses have already appeared elsewhere 



(Wilms et al. 19991; Wilms et al. 1999c). Results from the 



monitoring observations are presented in Paper II. 

The remainder of this paper is structured as follows. We 
start with a description of our RXTE data analysis method- 
ology (^). We then present the results from the spectral 
analysis from RXTE (§|) and ASCA (§§). The RXTE tim- 
ing analysis is discussed in section |H| We discuss our results 
in the context of current physical models for the soft state 
(§[]) and then summarize the paper (§[?]). 



2 RXTE DATA ANALYSIS METHODS 

Onboard RXTE are two pointed instruments, the Propor- 
tional Counter Array (PCA) and the High Energy X-ray 
Timing Experiment (HEXTE), as well as the All Sky Moni- 
tor (ASM). We used the standard RXTE data analysis soft- 
ware, ftools 4.2. Spectral m odeling was p erformed mostly 
with XSPEC, version lO.OOz flArnaud 199r| . We used essen- 
tially the same data extraction and analysis strategy as in 



our analyses of RXTE data from GX 339-4 and V1408 Aql 
(4U1957+11); therefore, we only present a brief summary of 
this strate gy. For detailed inform ation we refer to Wilms et 
al. ( |l999a| ) and Nowak & Wilms ( |l999| ). 

The PCA consists of five co-aligned Xenon (with an up- 
per Propane layer) proportional counter units (PCUs) with 
a total effective area of about 6500 cm 2 . The instrument is 
sensitive in the energy range from 2 keV to ~ 100 keV ( Ja 



hoda et al. 1996), although the response matrix is best cal- 
ibrated in the energy range of w 2. 5-30 keV. We chose to 
increase the signal to noise ratio of our data by analyzing 
top Xenon layer data only. Background subtraction of the 
PCA data w as performed analo gously to our previous study 
of Cyg X-l ( Dove et al. 1998 ). To reduce the uncertainty 
of the PCA background model, we ignored data measured 
in the 30 minutes after South Atlantic Anomaly (S AA) pas- 
sages. The XSPEC 'corrfile' facility was used to renormalize 
the background file for all observations, typically by 1-2%. 

For short (Js 10 ks) segments of observations, we as- 
sumed the uncertainty of the data to be purely from a Pois- 
son distribution. For analysis of the full (> 100 ks) observa- 
tions, we accounted for the remaining detector calibration 
uncertainty by using the en ergy de pendent systematic errors 
described by Wilms et al. (1999a). 

HEXTE consists of two clusters of four Nal/Csl- 
phoswich scintillation counters that are sensitive from 15 
to 250 keV. A full des criptio n of the instrument is given 
by Rothschild et al. ( 199£ ) . Background subtraction is 
done by source-background switching. We used the standard 
HEXTE response matrices of 1999 August, and considered 
data measured above 20 keV. 



3 RXTE SPECTRAL ANALYSIS 
3.1 LMC X-3 

RXTE observed LMC X-3 during three observing blocks on 
1996 November 30, 1996 December 2, and 1996 December 4 
and 5. A log of the observations can be found in Table |l[ 
During part of the observations, some detectors were turned 
off; however, we added all good data together and combined 
all background models. The total response matrix was the 
average of the respective estimates for each individual detec- 
tor combination, weighted by the fraction of photons mea- 
sured during the time that such a combination was active. 

The lightcurve of LMC X-3 during the long observation 
is displayed in Fig. During the first part of the observa- 
tion the source is characterized by a slight increase of the 
PCA count rate from ~ 80 to ~ 85 cps per PCU. The source 
was at a relatively steady level of ~ 77 cps per PCU during 
the last part of the observations. This variation was also ob- 
served in the HEXTE (Fig. [j]). To see whether this variation 
was also detectable spectroscopically, we first analyzed sep- 
arately the observations indicated in Table [j]. For this anal- 
ysis, we ignored data below PCA channel 7 and above PCA 
channel 50, which corresponds to an energy bandpass from 
3.6-20 keV. (Here and throughout our spectral discussions, 
PCA channel number refers to the 'standard_2f channels, 
1-129, as read by XSPEC.) 

We first modeled the spect rum via the standar d multi- 
temperature disc black body (Mitsuda et al. 1984) plus a 
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Table 1. Observing log for the long RXTE observation of LMC X-3. 



Obs. Date Exposure Count Rate 





JD-2450000 


ymd 


sec 


counts s 1 


a 


417.905 


1996.11.30, 09:43 


13800 


399.3 ±0.2 


b 


418.238 


1996.11.30, 17:42 


10200 


412.9 ±0.2 


c 


419.843 


1996.12.02, 08:14 


14000 


422.3 ±0.2 


d 


420.168 


1996.12.02, 16:02 


12700 


411.3 ±0.2 


e 


421.818 


1996.12.04, 07:38 


12400 


369.9 ± 0.2 


f 


422.100 


1996.12.04, 14:25 


11300 


388.8 ±0.2 


g 


422.374 


1996.12.04, 20:58 


5500 


384.2 ±0.3 


h 


422.660 


1996.12.05, 03:50 


16900 


374.2 ±0.2 


i 


422.960 


1996.12.05, 11:02 


11300 


334.9 ±0.2 


j 


423.238 


1996.12.05, 17:42 


16400 


395.0 ±0.2 


k 


423.502 


1996.12.06, 00:02 


15200 


392.2 ±0.2 



Exposure times shown are rounded to the closest 100 sec. The count rate is the total PCA background subtracted count rate. 



Table 2. Results of Spectral Fitting to the LMC X-3 Data. 



Obs. kT[ n ^disc ^compps T 

keV 10~ 2 



total 


1 


l/l+O 
24_ 


00 
01 


33.6I 09 






2 


n+0.1 
°-0.2 


total a 


1 


20±g 


01 
02 


37.3+J-* 






2 


g+0.2 
D -0.3 


a 


1 


231° 


00 
00 


25.61™ 






2 


q+0.0 
a -0.0 


b 


1 


24+°, 


00 
00 


28.2!°;° 






2 


q+0.0 


c 


1 


251°, 


00 
00 


2Q q+0- 6 






3 


Q + 0.0 


d 


1 


25+« 


00 
00 


9Q 7+0-6 
zy - '-0.5 






3 


q+0.1 
3 -0.1 


e 


1 


24 +0 


00 
00 


30.51°;° 






3 


4+0.1 

^-0.1 


f 


1 


24 +0 


00 
00 


30.91°;° 






3 


4+0.1 

4 -0.1 


g 


1 


24+° 


01 
01 


30.5l°;« 






3 


4+0.1 

4 -0.1 


h 


1 


231° 


00 
00 


30.8l°-j 






3 


q+0.1 
°-0.1 


i 


1 


231° 


00 
00 


30.8l°; 7 






3 


i+O.l 
i -0.1 


j 


1 


241° 


00 
00 


29.7l°j 






3 


i+O.O 
L -0.1 


k 


1 


241° 


00 
00 


28.8±g;| 






3 


,+0.0 


total b 


1 


161° 


04 
02 




2 


c+0.2 
°-l.l 






a 


1 


11+° 
L1 -o 


01 
01 




5 


g+0.4 






b 


1 


141° 


01 
01 




4 


g+0.6 
D -0.7 






c 


1 


161° 


01 

01 




2 


o+0.4 
8-0.2 






d 


1 


161° 


00 
00 




3 


Q + 0.1 

°-0.4 






e 


1 


171° 


00 
01 




1 


q+0. 3 






f 


1 


161° 


00 
00 




2 


Q+0.1 

a -0.5 






g 


1 


161° 


01 

01 




2 


1+0.3 
L -1.0 






h 


1 


151° 


01 
00 




2 


o+O.O 
8-0.4 






i 


1 


161° 


01 
01 




1 


c+0.0 
°-0.1 






j 


1 


161° 


00 
01 




2 


4+0.2 






k 


1 


151° 


01 
00 




3 


1+0.2 
L -0.1 











l/compps 


^Line 


EW 


X 2 /dof 


10 _1 


keV 


10~ 2 


10~ 4 


eV 




2 4+ 1 ' 1 








10 


50.4/47 


1 2+ 0,7 
1 - Z -0.6 






s-iili 


134 


37.8/46 


4 8+°' 4 






14 1+ 2 - 9 


22 


68.1/39 


q 7 +0.4 






1 6+°- 4 


24 


66.6/39 


3-8121 






2 0+°- 3 
z - u -0.3 


30 


99.6/39 


5-o±g:S 






1 S+0- 4 
1-8-0.3 


28 


61.8/39 


4 6+0-7 

4 - D -0.6 






1 9+°- 3 


32 


68.6/39 


44+0.7 
4 - 4 -0.6 






i r+0.3 
1 -°-0.4 


25 


59.3/39 


4 4+1-1 






1 6+°- 5 


28 


43.5/39 


4 0+°- 5 






1 4+0-3 


23 


128/39 


q q+0. 6 

•"-0.5 






1 O+0.3 
1 - a -0.4 


30 


54.5/39 


4 1+0-4 
^■ 1 -0.4 






1 q+0- 3 


31 


63.7/39 


4 8+°' 5 






1 4+0-3 
i - 4 -0.3 


23 


123/39 




621^ 


6 1+ ' 1 

°- -0.2 


10 1+ 3 ' 6 


170 


39.5/46 






9-5lE;:o 


13.81H 


230 


57.2/39 




48±* 


8.3l°; 2 


13.5+i- 7 


213 


46.4/39 




64l^ 4 


6 9+°- 2 


14.611;? 


228 


88.7/39 




43U 4 


5 q+0.1 
O.O_ 2 


13.711;° 


223 


71.7/39 




84±f 5 


45+0.1 

*-°-0.2 


13 1+ 1 ' 5 


228 


80.8/39 




40+1 1 


4 q+0-1 


11.0+ 2 ° 


190 


67.9/39 




50l?° 


4 2+0.2 


10.811;° 


190 


57.3/39 




4411 1 


4 6+ al 


ii.4ii;i 


202 


113/39 




89l?° 


5 6+°- 2 


12.2+1-, 


213 


68.8/39 




64l« 


6 0+ ' 1 


12.411'j 


209 


80.1/39 




sol? 


°-°-o.i 


1 2 2+ 1 ' 3 


210 


130/39 



Tin, -4disc : Peak multi-temperature disc temperature and normalization. T: Photon index of the power law. Apl: Power law 
normalization (photons kcV~i enr 2 s _1 at 1 keV). Acompps: Normalization of compps model. kT e : Coronal electron temperature for the 
compps model. t; com pp S : Compton y-paramcter for the compps model. ALine: Line normalization (photons cm -2 s^ 1 in the line). The 
Gaussian line was fixed at 6.4 keV with a width of a = 0.1 keV for the multi-tcmpcraturc disc plus black body models, or a = 1.4 keV 
for the compps models. EW: line equivalent width. Uncertainties are at the 90% confidence level for one interesting parameter 
(Ax 2 = 2.71). The interstellar equivalent column was fixed at Nn = 3.2 X 10 20 cm -2 . a Paramctcrs for a broad iron line with 

a = 1.4l°' 2 keV. b Parameters for a broad line with a = 1.4l°' 2 . 
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JD-2450000 

420 421 422 



o 



< 
o 

0- 



X 
m 
X 
H 

m 



1 2 3 4 5 6 

1996 December 

Figure 1. Background-subtracted RXTE PC A lightcurve of 
LMC X-3, rebinned to a resolution of 64 s (left y-axis), and to- 
tal background subtracted HEXTE count rate (sum of HEXTE 
clusters A and B; rightw-axis). The dashes indicate the individual 
observations of Table 111 



power- law component. We modeled absorption in the in- 
tervening interstellar medium by fixing the equivalent hy- 
drogen column to the value found from radio observations 
7V H = 3.2 x 10 20 cm" 2 (Staveley-Smith 1999, priv. comm.), 
and by usi ng the cross sections of Baluciriska-Church & Mc- 
Cammon (1992). Our preliminary fits to this model had 
strong residuals in the region around ~ 6.5 keV (Fig. ^| 
middle). We therefore included a narrow emission line at 
6.4 keV. We will return to the issue of the reality of this line 
feature below. Table ^ lists the parameters of the best fit 
models. 

The characteristic temperature of the phenomenologi- 
cal accretion disc component, fcTj n , can be determined to 
high accuracy (see also Paper II). The formal uncertainty of 
the values for kT- m presented in Table |^ are typically smaller 
than 0.005 keV. This rather small formal uncertainty value 
might indicate that the systematic uncertainty of the re- 
sponse matrix needs to be taken into account even for these 
short observations. These systematic uncertainties, however, 
are of roughly the same magnitude as the Poisson errors 
used in our analysis. This fact makes a formal error analysis 
rather difficult, and it is not clear whet her th e parameter 
error analysis methods of Lampton et al. ( 197C ) yield mean- 
ingful results in such a case. 

During the week over which the observations were per- 
formed, kT ln remained more stable than the other fit param- 
eters, although these other parameters had relatively weak 
variations. It has been observed that on much longer time 
scales (months to years) variations of LMC X-3 are corre- 
lated predominantly with strong variations of fcTj n , while 



0.10 
0.08 

I 

J: 0.06 

0.04 

0.10 
0.08 

I 

J 0.06 
0.04 



- 1.18 

^ :i.i6£ 

i : J3 

T - 1.14 w 

: I 

H 

- 1.12 * 
1 1.10 



418 419 420 421 422 423 424 
JD-2450000 



1.10 



1.12 1.14 

kT dist (keV) 



1.16 



Figure 2. Top: Peak disc temperature and Compton y-parameter 
as a function time for the observations of LMC X-3 presented in 
Tab. ^. Bottom: Peak disc temperature vs. Compton y-parameter. 



the phenomenological disc normalisation, Adi sc , remains rel- 
atively constant (Paper II, Ebisawa et al. 1993). In the lit- 
eral interpretation of the multi-temperature disc black body, 
Adisc oc r 2 n cosi, where n n is the inner radius of the accre- 
tion disc and i is its inclination. The stability of Adi sc thus 
has been interpreted as evidence that the inner radius of the 
accretion disc stays remarkably constant on long tim escale s 
( Ebisawa et al. 1993 ). As discussed by Merloni et al. ( 1999 ), 
however, this is not necessarily true for a more physical disc 
model that accounts for radiative transfer effects, Doppler 
blurring, and gravitational redshifting. For example, a model 
wherein only the accretion rate and not the inner disc ra- 
dius varies formally can be fit with a disc black body with 
varying yldisc (Merloni et al. 1999). 

As a more physically motivated description, we have 
also fit the Compton scattering model of Poutanen & Svens- 
son (1996) (compps). Specifically, we simulate a disc black 
body Compton upscattered via a slab geometry corona with 
fixed unity covering fraction over the disc. A purely ther- 
mal coronal electron distribution was assumed, and the 
width of the Fe line was fixed to be a ~ 1.4 keV (see be- 
low). We then fit the peak temperature of the disc black 
body, as well as the electron temperature and Compton y- 
parameter of the corona. Results are presented in Table [j] 
and in Fig. ^| Here the peak temperature of the disc is seen 
to vary more strongly over the course of the observation 
than for the multi-temperature black body plus power-law 
models. The variations in the lightcurve are consistent with 
being predominantly driven by variations of the Compton 
2/-parameter, with the large count rate change from observa- 
tions a to e corresponding to a drop in y. The slight increase 
in count rate seen in observations h-k corresponds to a slight 
increase in y. There is a tendency for the peak disc tem- 
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100.0 E 



10.0 = 



1.0 - 



0.1 - 



2 
1 

Otr 
-1 
-2 



2 - 
1 

-1 

-21 



10 

Energy (keV) 



20 



30 



Figure 3. Top: Observed spectrum and disc black body plus 
power law plus broad iron line model fit to the long observation 
of LMC X-3. Dashed lines show the individual model components 
folded through the detector response. Middle: Residuals in units 
of a between the data and the model, with a narrow iron line. 
Bottom: Residuals for the best fit model, including a very broad 
iron line. (Residuals for the compps model are virtually identical.) 



perature and Compton {/-parameter to be anti-correlated; 
however, this anti-correlation is dominated by observations 
a and b. 

Note that the lightcurves cover approximately 3, 1.7 day 
orbital periods of LMC X-3. Recent work with these (now 
archived) RXTE observations, along with more recent long- 
term pointings, reveal evidence for a weak (few percent) 
orbital modulation of the X-ray source if one folds the 
lightcurves on the known orbital period (Boyd & Smale, 
in prep.). Such variations would not be unexpected as pre- 
vious observations have suggested a near edge-on s ystem to 
accoun t for the ellipticity of the optical lightcurve (van dei 
Klis et 



al. 1983). The nature of these variations, possibly 



due to weak scattering, will be discussed further by Boyd & 
Smale (2000). 

We now consider fits to the summed spectrum from 
the long observation. Taking PCA channels 7 to 58 (~ 3.6- 
26keV), we again fit two models: a disc black body plus 
power law plus gaussian line, and a compps model plus gaus- 
sian line. The best fit disc black body plus power law models 
are shown in Fig. |§| For both classes of models, we formally 
required a broad Fe line component with equivalent width 
> 150 eV with widths of a » 1.4 keV (Table |). This is ap- 
proximately twice the PCA spectral resolution. Despite the 
low reduced \ 2 1 we do not consider such a strong line to 
have been detected definitively in LMC X-3 as there is some 
concern about systematic errors in both the PCA response 
and the spectral model. 

Excluding the line, residuals in the line region are ~ 1%, 
similar to the systematic uncertainties in the PCA response. 
In addition, for the models that we have fit, energy channels 



6 keV are dominated by the disc component, whereas ener- 
gies "it lOkeV are dominated by the power-law/Comptonised 
component. The crossover region corresponds to the Fe line 
region (see Fig. |^). Thus the line strength is especially sensi- 
tive to any errors in modelling the continuum spectral shape, 
such as, e.g., the approximation made by the diskbb+power 
model, where the comptonised spectrum is represented by 
a pure power-law even at energies where the comptonised 
photons have energies comparably to the seed photon en- 
ergy. This is in contrast to AGN, for example, where the 
continuum through the line/edge region is thought to be 
reasonably well-approximated by a featureless power law, 
and therefore is not subject to systematic uncertainties in 
the continuum model. 

For the case of AGN, the reality of the line can also be 
verified (independently of response matrix uncertainties) by 
taking the ratio of the (nearly power-law) observed spectrum 
to that of the power-law spectrum of the Crab nebula. In 
principle, a similar procedure could also be used in case of 
LMC X-3. Simulating a compps spectrum both with and 
without a line, utilizing the fit parameters for LMC X-3 
discussed above, shows that a division of these two spectra 
would reveal the presence of a line. There exists no template 
spectrum, however, with which to do this comparison. We 
caution, therefore, that the LMC X-3 line in reality may be 
weaker and/or narrower than that obtained here. 

The Comptonisation model implies a coronal temper- 
ature of w 60keV, and a Compton {/-parameter of y = 
AkT CB /mc 2 r os » 0.06, which implies an optical depth of 
t os ~ 0.1. Due to the limited energy range of our obser- 
vations, however, these fits are not unique. We assumed 
a purely thermal electron distribution, but we note that a 
purely non-thermal Comptonizing electron distribution fits 
the data equally well. Data at energies significantly higher 
than 50keV are required to break this degeneracy. 

As regards the power-law component, when consider- 
ing the HEXTE data as well, we detect the source out to 
ps 50keV, the highest energy at which this source has ever 
been detected. The HEXTE data, however, do not place 
much stronger constraints than were obtained with the PCA 
data alone. Including the HEXTE data, the observed power 
law is r = 2.81q'2- We do not detect any noticeable curva- 
ture in the spectrum, but our limits on the presence of such 
curvature are weak. 



3.2 LMC X-l 

Historically, LMC X-l has shown relatively little long time 
scale variability (Syunyaev et al. 1990, see also Paper II). 
Consistent with this, LMC X-l showed no obvious variations 
on day long time scales during the long RXTE observation 
of 1996 December 6 to 8. The combined spectrum is a typical 
example of the spectrum of LMC X-l (see Paper II). Taking 
our screening criteria into account, we obtained 80ksec of 
data with 5 PCU on, and another 30ksec of data wherein 
there were 3 or 4 PCU on. Again, these data have been 
added together as for LMC X-3. This resulted in a detection 
of LMC X-l out to 20 keV in the PCA. In HEXTE, LMC X- 
1 is detectable out to the same energy, but due to its low 
count rate (the total HEXTE count rate is 1 counts s _1 ), 
no spectral information can be obtained. 

We modeled the data (PCA channels 7 to 52, » 3.6- 
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Table 3. Results of Spectral Fitting to the LMC X-l Data. 

kT in A disc Acompps r A PL kT a Scompps a A Linc EW x 2 /dof 

keV 10~ 2 1CT 1 keV 10" 2 

0-9llooi 58.0+" ... 3.1+'" 2.1+°-^ 0.1 7±* 38 71.8/41 



-0.01 ""'"-3.0 ■•■ "• J --0.1 0.6 

0.88±°-°i 76.3±|S ... 2.9±oi L5 +o.4 0.98+";^ 32+J" 195 46.4/40 



-0.01 lu - J -6.0 ■•■ "-0.1 x - J -0.5 

°- 86 ^ao2 •■• L9 to.3 ■•■ •■■ 68 -l 4 - 4 -o'.3 °- 97 ±o:i4 33 -l° 207 48.5/40 
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Figure 4. Same as Fig. 3 for the long observation of LMC X-l. 



22 keV) with the same spectral models as that for LMC X- 
3, using an equivalent hydrogen column of Nu = 7.2 x 
10 21 cm -2 (Staveley-Smith 1999, priv. comm.), taken from 
21 cm measurements. Again, the data required the presence 
of an iron line at 6.4 keV. The resulting best fit parameters, 
with a narrow and with a broad iron line, are shown in Ta- 
ble [j. In Figure ^, we display the residuals of the disc black 
body fits, showing the improvement of assuming a broad 
iron line. As for LMC X-3, without a broad line the rela- 
tive deviation between the data and the model in the line 
region is ~ 1%. Again, the line region corresponds to the 
transition region between the disc black body and power- 
law/Comptonized emission components of the spectrum, so 
caution is warranted in interpretting the fit parameters of 
the line. 

Our spectral parameter values are in general agreemen t 



with earlier data (Ebisawa et al. 1989; Schlegel et al. 1994) 



although our power-law index appears to be rather soft com- 
pared to the earlier measurements. As we show in Fig. ^, the 
soft (r ~ 3) power-law component dominates the observed 
flux above ~7keV, but also contributes an appreciable frac- 
tion of the total soft flux. Thus the lower energy channels 
are also in part determining the value of the fitted spec- 
tral index. In terms of the compps Comptonisation model, 
if we choose a purely thermal, slab-geometry corona with 
unity covering factor, we find similar coronal parameters 



as for LMC X-3; however, the disc temperature is some- 
what lower so the observed spectrum is more dominated by 
Comptonized photons than that for LMC X-3. Again, purely 
non-thermal models are also permitted, indicating that the 
coronal geometry and parameters are not uniquely deter- 
mined by these data. 

Since our observations cover a larger energy range than 
many previous missions, they ensure that our best fit models 
at least describe the s pectr al shape quite well. In accordance 
with Schlegel et al. (1994), we do not find evidence for a 



broad iron edge feature in the data. Such a feature might be 
artificially introduced in the data for instruments where the 
upper energy boundary is so low that the power-law index 
cannot be well constrained. 



4 ARCHIVAL ASCA DATA 

The High Energy Astrophysics Archive (HEASARC) con- 
tains one ASCA observation of LMC X-l (1995 April 2) 
and two observations of LMC X-3 (1993 September 23, 1995 
April 15). To the best of our knowledge, these observations 
have not been published previously. We extracted the data 
from all four instruments on ASCA, the two solid state de- 
tectors (SIS0 and SISl) and the two CIS detectors (GIS2 
and GIS3); however, here we only discuss the SIS data, as 
the SIS has a more reliable low energy response and a better 
spectral resolution. 



We used the SISCLEAN tools ( pay et al. 1998; ), with 
the default values, to remove hot and flickering pixels. Fur- 
thermore, we filtered the data with the same cleaning crite- 
ria outlined by Brandt et al. (1996); however, we took the 
more conservative values of 10° for the minimum elevation 
angle and 7 GeV/c for the rigidity. We obtained background 
estimates by extracting source-free regions near the edges of 
the SIS chips. The spectra were rebinned so that each en- 
ergy bin contained a minimum of 20 photons, and we only 
fit SIS data in the 0.5 to lOkeV range. Cross-calibration 
uncertainties between the two SIS detectors were accounted 
for by introducing multiplicative constants (always found to 
be within ^ 1% of each other) for each detector in all fits. 

N eutral hydrogen abso rption of a disc black body spec- 
trum (Mitsuda et al. 1984) provides a rough description of 



the ASCA data. Such models yielded reduced x 2 values 2 
for LMC X-l and w 1 for LMC X-3. All three observations, 
however, exhibited evidence for excess emission at energies 
<; 5keV. Adding a power law (photon flux oc E~ r , with 
r ~ 2.2-2.5) improved the reduced x 2 to ~ 1.1 and ~ 0.7 
for LMC X-l and LMC X-3, respectively. As for the RXTE 
observations, we also were able to fit the compps model; how- 
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Table 4. Results of fitting the LMC X-l and LMC X-3 ASCA data. 
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See Table. ^ for an explanation of the symbols. fcT ray , A ra y are the temperature and normalization of the Raymond-Smith plasma 



component (Raymond 



Smithl977). Parameters typeset in italics were frozen at the indicated value. Errors are at the 90% confidence 
level for one interesting parameter (Ax 2 = 2.71). 
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Figure 5. ASCA observations of LMC X-l. Top panel shows 
the compps plus a Raymond-Smith plasma model fit to the data, 
with N H fixed to 7.2 x 10 2 1 cm . Dashed lines show the indi- 
vidual model components folded through the response matrix of 
the SIS0 detector. Lower panels, from top to bottom, show data 
residuals for: compps model with fixed Nn and no plasma compo- 
nent, compps model with fixed iVjj and a Raymond-Smith plasma 
component, and compps model with free Nn and no plasma com- 
ponent. 



ever, given that the ASCA data cuts-off at 10 keV there was 



not enough leverage to fit both kT c and y c , 



We there- 



fore froze the coronal electron temperature to 50 keV for all 
three observations. Results for these fits are presented in 
Table g and Fig. |. 

Both LMC X-3 observations are at very low flux values. 
We compared the 3-9 keV flux values from these observations 
to those obtained from our RXTE monitoring observations 
in Paper II. Using an SIS/PCA normalisation ratio of 0.7, 
we estimate that the ASCA observations would correspond 



to ASM count rates of ~ 1.3 cps and ~ 0.6 cps. As discussed 
in Paper II, these low count rates occur near transitions from 
the soft to hard state of LMC X-3. We note that the lower 
flux ASCA observation of LMC X-3 shows a significant drop 
in peak disc temperature, consistent with the trends seen in 
the state transitions discussed in Paper II. 

Our main purpose for examining the ASCA data is to 
determine if there is evidence for an Fe line in the spec- 
trum and an Fe-L complex near energies of m 1 keV. Line 
features near 1 keV are a common occurrence in photoion- 
ized plasmas close to sources emitting hard X-rays (e.g., in 
eclipse in Vela X-l, Nagase et al. 1994). We would also ex- 
pect such features in models of warp ed ac cretion discs with 
winds, similar to those of Schandl ( 1996[ ), or in models of 
wind-driven limit cycles (Shields et al. 1986), as might be 
relevant for producing the long term periodicity of LMC X-3 
(Paper II). Such a line complex also was seen in th e spec- 



trally very sim ilar source V1408 Aql=4U1957+ll (Nowak 



Starting with the compps models, none of the three 
ASCA observations show evidence of an Fe line. If we freeze 
the line energy at 6.4 keV or 6.6 keV and the line width at 
0.3 keV, we find an upper (90% confidence) limit of 60 eV for 
the line equivalent width for the LMC X-l observation. Nar- 
rower lines (down to the ASCA resolution of ~ 0.1 keV) pro- 
duce more stringent equivalent width limits. For the slightly 
fainter 1995 April 15 and 1993 September 23 LMC X-3 ob- 
servations, the equivalent width 90% confidence limits are 
90 eV and 320 eV, respectively. 

The reduced % 2 values with the compps model are sig- 
nificantly less than one for both LMC X-3 observations, and 
therefore we do not attempt to fit a plasma component to 
these data. The compps model fit to the LMC X-l observa- 
tion yields a reduced x 2 ~ 1; however, this fit requires an 
Nn value 12% lower than the estimates of Staveley-Smith 
(1999, priv. comm.). Fixing the Nn column density to this 
value worsens the reduced x 2 to ~ 2 and yields significant 
residuals in the Fe-L complex region (Fig. ^). This residual 
can be removed and the reduced x 2 value returned to ~ 1 
by adding a plasma component with kT ~ 0.8 keV (here 
modeled with the XSPEC raymond model, after Raymond 
& Smith 1977). Thus the reality of any Fe-L complex is seen 
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Figure 6. Fourier frequency times Power Spectral Density (PSD) 
for LMC X-l in the 0-3.3 keV (solid diamonds, lowered by a 
decade), 3.3— 4.7keV (clear diamonds), 4.7-9.1 keV (solid squares, 
raised by a decade) energy bands. PSD normalisation is such that 
integrating over positive frequencies yields the mean square vari- 
ability divided by the square of the mean for the lightcurve an- 
alyzed. Lines represent the expected level of positive 1-<t noise 
residuals after subtracting Poisson noise from the PSDs. 



to depend strongly upon accurate modeling of the neutral 
hydrogen absorption of this source. Differing metal abun- 
dances in the LMC also could lead to errors in interstellar 
absorption models, adding further uncertainties to the level 
of any plasma component. If real, however, a plasma com- 
ponent at the levels allowed by the ASCA data with a fixed 
neutral hydrogen column would be readily detectable by the 
X-ray Multiple Mirror Mission, which would also resolve the 
complex into individual line components (Note also that the 
residuals seen at ~ 1.7 and ~ 2.2 keV are likely related to 
well-known systematic features in the SIS responses matri- 
ces.) 



5 TIMING ANALYSIS 

We employed Fourier techniques, in t he same manner as fo r 
our RXTE observations of Cyg X-l flNowak et al. 1999aj), 



to study the short timescale variability of LMC X-l and 
LMC X-3. Specifically, we used the same techniques for es- 
timating dea dtime corrections ( Zhang et al. 1995| ; Zhang & 
Jahoda 1996) to the Power Spectral Density (PSD), and for 
estim ating uncertainties a nd the Poisson noi se levels of the 
PSD (iLeahy et al. 19831 Ian der Klis 19891). In the anal- 



ysis discussed below, we use lightcurves with 2~ s resolu- 
tion constructed from the PCA top Xenon layer data only. 
Furthermore, we subdivide the lightcurves into three energy 
channels: 0-3.3 keV, 3.3-4.7 keV, and 4.7-9.1 keV (absolute 
PCA channels 0-8, 9-12, 13-24, respectively). We chose 
these energy ranges as they are relatively background free, 
and furthermore the three channels have roughly equal count 



rates. The lowest energy channel is dominated by the disc 
component, whereas the highest energy channel is sampling 
the power law component (iQ). 

To maximize the signal to noise ratio, we created a 
single PSD in each energy channel that was averaged over 
the entire duration of our observations. We further aver- 
aged PSDs constructed from data segments of 1024 sec du- 
ration (76 segments for LMC X-l, 65 segments for LMC X- 
3), and we logarithmically binned the PSD over frequencies 
/ — > 1.3/. For the case of LMC X-3, weak variability with 
root mean square (rms) amplitude of 0.8% of the mean was 
detected in the 10 -3 -10 -2 Hz range. This is consistent with 
background fluctuations. 

Significant variability above the noise, however, was de- 
tected in LMC X-l. Root mean square variabilities of 6.5%, 
6.3%, and 7.1%, lowest to highest energy channel, were 
found for the frequency range 7 x 10 _4 -0.3 Hz. The PSDs for 
these energy channels are shown in Fig. [| All three PSDs 
are approximately proportional to / _1 from 7x 10 _4 -0.2 Hz. 
The highest energy band PSD, however, shows evidence for 
a rollover at ~ 0.2 Hz. Note that for this latter PSD there are 
only three or four points at frequencies higher than 0.2 Hz 
that sharply decline, followed by several points with a more 
gradual decline roughly proportional to / _0 ' 5 . (Fig. ^ shows 
the PSD multiplied by frequency.) The f~ ' 5 behaviour is 
also readily apparent at high frequencies in the PSD of the 
low energy channel as well as possibly in the PSD of the 
middle energy channel. An / _0 ' 5 proportionality is exactly 
that expected for positive noise residuals if the mean Pois- 
son noise level was slightly underestimated (see Nowak et al. 
1999a). The uncertainty in the noise level makes it difficult 
to assess the frequency-dependence of the rollover as a func- 
tion of energy. The 4.7-9.1 keV energy band PSD values at 
frequencies 0.2-2 Hz, however, clearly lie below an extrap- 
olation of the low-frequency PSD behaviour. Removing the 
break altogether would require that we overestimated the 
Poisson noise level, which seems doubtful given the oc f~ ' 5 
behaviour of the noise-subtracted PSD residuals. 

These results for LMC X-l are roughly consistent with 
the previ ous RXTE observations reported by Schmidtke et 
al. (1999) in terms of ove rall P SD amplitude and shape. 
However, Schmidtke et al. (1999) did not attempt any noise 
subtraction and did not search for breaks in any of the PSD. 
We do not find any evidence for a 0.08 H z qua si-periodic 
oscillation as reported by Ebisawa et al. (1989 ), wh ich is 
consistent with the results of Schmidtke et al. ( 1999). The 
0.08 Hz QPO reported by Ebisawa et al. (1989), however, 



is nearly coincident with the expected level of the residuals 
after noise subtraction. Thus the previously reported QPO 
may have be an artifact of a misestimation of the Poisson 
noise level, as opposed to the lack o f detection here and in 
the work of Schmidtke et al. ([l99§) being due to an inter- 
mitant nature of such low- frequency QPO. 



Contrary to the results reported by Treves et al. (1988) 
we do not detect any variability from LMC X-3. We note 
however, that the PSD reported by Treves et al 
the characteristic f~ 
noise residuals. We 



1988) has 



shape and amplitude expected from 
believe, therefore, that the results of 
Treves et al. (1988) are consistent with those reported here. 

We further studied the temporal behaviour of LMC X- 
1 by computing the Fourier frequency-dependent time lags 
and variability coherence between the various energy bands 
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Figure 7. Coherence between the variability in the 0-3.3keV and 3.3-4.7keV energy bands (left), the 3.3-4.7keV and 4.7-9.1keV energy 
bands (middle), and the 0— 3.3keV and 4.7— 9.1keV energy bands (right). 



(see Vaughan & Nowak 1997, and references therein). Due 
to the low count rates, only weak l-cr upper limits of 
<; 2 sec (//10 -2 Hz) -1 could be placed on the time lags 
between the lowest and highest energy channels over the 
frequency range 10 _3 -0.3 Hz. We were able to measure the 
coherence among the variability in the various energy bands. 



As we 'lave previously discussed (Vaughan 
see also Bendat & Piersol 
7 2 (/), which is always < 



Nowak 1997; 
1986), the coherence function, 
1, is a measure of the degree of 



i ng partly via a wind with velocity v v <; 600- 1100 km s" 1 
( Hutchings et al. 1982 ; Hutchings et al. 1987 ) . This has a 
number of implications for the accretion flow. 

Taking a p resumed mass of 6 Mn and a wind velocity 
of 600 km s _1 flHutchings et al. 1987), the Hoyle-Lyttleton 
accre tion radius (Hoyle & Lyttleton 1939; Bondi & Hoyk 
1943) is then 



linear correlation between two time series, or alternatively, 
it is a measure of how well one time series can be predicted 
from the other. 

We present the coherence function for the LMC X-l 
variability lightcurve, logarithmically binned over frequen- 
cies / — > 2/, in Fig. u\. The lowest and middle energy 
band lightcurves are fairly well-correlated and have a coher- 
ence function between them of 7 2 (/) 0.9. Both of these 
lightcurves, however, are very incoherent with the highest 
energy band lightcurve, with 7 2 (/) ~ at 10 -3 Hz and then 
rising approximately as oc log(/). Due to uncertainties in 
the Poisson noise level there are additional systematic un- 
certainties above w 0.3 Hz; however, the coherence is clearly 
below the near unity levels that we typically observe during 
t he low/hard states of BHC such as Cyg X-l and GX339-4 
(|Nowak et al. 1999a|; |Nowak et al. 1999b|). 



2GM 



■ 4xlO n cm 



M 
6M F 



600 kms- 



which should be compared to the binary separation of 
a ~ 2.2 x 10 12 cm given a 4.2 d ay period and a 20 Mq 
companion ( Hutchings et al. 1987 ). The "circularisation ra- 
dius", however, will be much smaller, by a factor of the 
order (w or b/^w) 2 , where u or b ~ 260 km s _1 is the orbital 
velocity. This estimate comes from assuming the accreted 
angular momentum per unit mass goes as J?hl ■ v OT b, and 
then setting this value equal to the Keplerian value at the 
circularisation radius. Thus, we expect a disc circularisation 
radius of 8 x 10 10 cm (i; w /600 km s -1 ) -4 . The disc is thus 
expected to be smaller than would be usual for Roche lobe 
overflow, perhaps substantially so. For w w = 2000 km s _1 , 
the expected circularisation radius is ~ 750 GM/c 2 . 



& Taam 1985 ; 


Taam h 


z Fryxell 1988; 


Benensohn et al. 1997 


Ruffert 19971 



6 DISCUSSION 

6.1 Variability of LMC X-l 



foam fc Fryxe 



1 198c 



tuffcrt 199S| ) showing 



that wind dominated accretion inherently exhibits fluctu- 
ations in the net accretion rate onto th e compact object 



More rec ent 3-D numerical integrations ( fluffcrt 1997 



The mc st obvious striking difference between LMC X-l and 
LMC X-3 is that the former shows moderate (7% rms) vari- 
ability whereas the latter shows essentially none (< 0.8% 
rms) on timescales shorter than Iks, despite the fact that 
the X-ray spectra of these two objects are remarkably similar 
in terms of shape and overall flux. We hypothesize that the 
differences lie within their respec tive modes of accreti on. A s 
discussed by Cowley et al. (1983), van der Klis et al. (1983), 



Ruf- 



fert 199£) tend t o show smaller amplitu de fluctuations than 



and Cowley et al. (1994), LMC X-3 is consistent with accret- 
ing via Roche lobe overflow. One would therefore expect the 
accretion rate of LMC X-3 to be relatively steady (although 
see the discussion of the long-term periodicity presented in 
Paper II). LMC X-l, on the other hand, is most likely accret- 



earlier 2-D work (Benensohn et al. 1997); however, accretion 
rate fluctuations of O(10%), in rough agreement with the ob- 
served variability, are expected. LMC X-3, being Roche lobe 
fed, would not be expected to exhibit these fluctuations. 

A disc with a truncated outer edge may also yield an 
explanation of the possible rollover in the high-energy band 
PSD seen at ~ 0.3 Hz. For a = 0.3, H/R ~ L/L md ~ 0.3, 
M — 6 Mq, the viscous timescale at R — 750 GM/c 2 is ~ 
20 sec, which is in rough agreement with the time scale of the 
possible rollover. This is consistent with the idea that any 
accretion rate fluctuations on time scales shorter than the 
viscous time scale at the circularisation radius are smoothed 
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the 'power law component' (oc 3.6 1 7 /[l — 7]) in an anal- 
ogous way. The fastest variation occurs on time scales of 
« 0.4 days. Given a = 0.3, H/R ~ L/LEdd ~ 0.3, and 
M — 9 Mq, this corresponds to a viscous time scale at a 
for the low/hard states of Cyg X-l and GX 339—4 (Nowak radius of ~ 8 x 10 10 cm, approximately one quarter of the 



out as they propagate into the inner X-ray emitting regions 
of the disc. 

As discussed above, the low variability coherence values 
for LMC X-l at / > 10~ 3 Hz are unlike what we observed 



et al. 



Nowak et al. 1999b). We are aware of three pos- 



Nowak 



sibl e wa ys of explaining such low values (Vaughan , 
1997 ]) . ^ irst, the low and high energy bands may be com- 
pletely unrelated to each other. This seems unlikely as all 
three energy bands exhibit a roughly cx / - PSD with com- 
parable amplitudes. Second, the upper energy band PSD 
may be related to the lower energy band PSD via a non- 
linear transfer function. This is a viable option, given the 
complicated physics of "coronal formation" and the fact that 
non-linear transfer functions couple power from f requencies 
/ in o ne ban d to multiples of / in other bands (Bendat & 
Piersol 1986). This latter case is especially difficult to detect 
in the /~ PSD because of its scale-free nature. 

The third possibility is perhaps the most viable option. 
"Mixing" of independent components with widely varying 
intrinsic time lags can lead to the loss of coherence (see the 
discussion presented by Nowak et al. 1999b). For example, 
if one component which shows a correlation between hard 
and soft X-rays (e.g., local dynamical time scale disc fluc- 
tuations leading to and increase in both seed photons and, 
very slightly delayed, upscattered hard photons) is mixed 
with another, nearly equal in power but independent compo- 
nent with an anti-correlation between soft and hard X-rays 
(e.g., matter being transferred between disc and corona on 
thermal or viscous time scales, with comparably long delays 
between soft and hard photons), then a strong loss of co- 
herence is expected. In a simple two-component model the 
coherence function, 7 2 (/), is given by 

R 2 (f)+2R{f) cos(9 b -9 a ) 



7 2 (/) 



[R(f) + 1] 2 



(2) 



Here R{f) is the ratio of Fourier frequency-dependent cross 
correlation amplitudes for the two independent processes be- 
ing considered (soft vs. hard X-ray variability, process b over 
process a), and 9 a , 8 b are the Fourier frequency-dependent 
phase delays between the hard and soft X-ray v ariability for 
the ind ependent processes a and b, respectively ( Vaughan & 
Nowak |l997| ; |Nowak et al. 1999b| ). (Time delay, r = 9/2nf.) 
For the plots in Fig. @, if R(f) « 1 and 9 b - 9 a ~ tv at 
/ = 10~ 3 Hz, then 7 2 (/) ~ there as well. The increase 
in coherence towards higher frequency could be due to a 
decrease of either 9b — 9 a or R(f). The latter might be rea- 
sonably expected, for instance, if one process is happening 
predominantly on 0(1000 sec) viscous time scales (again, 
perhaps matter/energy transfer between disc and corona), 
whereas the other one is happening over a broader range of 
time scales (fluctuations in the seed/Comptonized photons). 



6.2 Spectroscopy of the Soft State 

The spectroscopy of LMC X-3 perhaps bears out the no- 
tion of there being an anti-correlation between the disc 
and coronal components on long time scales, as shown in 
Fig. ^. A similar result can be obtained with the more phe- 
nomenological disc black body plus power law models pre- 
sented in Table 2. The photon flux in the 'disc component' 
(oc AdiBckT^ lBC ) is anti-correlated with the photon flux in 



disc circularisation radius. It is therefore more likely that 
this variation occurs on a thermal or dynamical time scale 
closer towards the center of the system. Thermal time scale 
transfer of energy between disc and corona is again consis- 
tent with the results shown in Fig. [| 

Unfortunately, it is difficult to say much more about 
the spectrum other than it is well-fit by a disc black body 
plus power law model. Given a 70° inclination, 9 Mq black 
hole, and 50kpc distance for LMC X-3, and a 30° inclina- 
tion, 6 Mq black hole, and 50kpc distance for LMC X-l, the 
disc normalisations in Tables ti and H are roughly consistent 
with inner disc radii of ~ 6 GM/c . This is quite remark- 
able, given the simplicity of the disc black body model (cf., 
Merloni et al. 1999). As regards the power law tail, other 
soft state black hole candidates exhibit evidence for non- 



thermal Comptonisation in this tail (Glierliriski et al. 1999). 
Here, however, as we were only able to go out to 50keV and 
20keV for LMC X-3 and LMC X-l, respectively, we were 
able to describe the data completely adequately with pure 
thermal Comptonisation. 

The extent to which a strong, broad Fe line is required is 
likewise still somewhat uncertain. The w 150 eV equivalent 
width, a ps 1.4 keV line is consistent with expectations from 
the basic picture outlined above: a disc extending down to- 
wards the black hole horizon, with a large covering fraction 
corona sitting on top of it. (A slight Fe overabundance, or 
perhaps beaming of the coronal radiation towards the disc 
might also be required for the implied equivalent widths.) 
However, as the Fe line occurs in a region where the two basic 
spectral components — soft disc spectrum and hard power 
law/Comptonisation spectrum — cross one another, the Fe 
line parameters are especially sensitive to our assumptions 
regarding the underlying continuum model. Again, this is in 
contrast to AGN where it is presumed that the underlying 
continuum in the Fe line/edge region is well-approximated 
by a simple power law. 



7 SUMMARY 

In this work, we have presented analysis of long (170ksec) 
RXTE observations of LMC X-l and LMC X-3, as well as 
analysis of archival ASCA observations of these objects. The 
primary results of these analyses are as follows. 

• Both LMC X-l and LMC X-3 are well-fit by a simple 
spectroscopic model consisting of a disc black body and a 
power law. The power law, which was seen out to 50keV 
for LMC X-3, shows no obvious curvature. The data are 
consistent with exhibiting a broad Fe Ka line with equiv- 
alent width ;$ 150 eV; however, the line could in reality be 
narrower and weaker. ASCA observations do not reveal the 
presence of a narrow Fe Ka line, and the presence of any 
Fe L-complex in LMC X-l is almost completely dependent 
upon the modeling of the neutral hydrogen column towards 
that source. 

• As a more physical interpretation of these results, the 
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spectra are consistent with Comptonisation of a disc spec- 
trum by a corona with electron temperature kT c ~ 60 keV 
and Tes ~ 0.1. These fits are not unique, however, as the data 
do not extend beyond 50 keV. Purely non-thermal Compton- 
isation models are also permitted. 

• LMC X-3 shows no variability on 1 ks time scales; 
however, it does show variability on half day time scales and 
longer. The time scale of these variations are slightly faster 
than the viscous time scales of the outer disc. Spectral mod- 
els of these variations are consistent with an anti-correlation 
between the disc and Compton corona components of the 
spectrum 

• LMC X-l shows rms variability of a few percent on 
faster than 1000 s timescales. This variability, along with 
with a possible power spectrum rollover on time scales ^ 5 s, 
is consistent with accretion rate variations in a wind fed 
system where the disc is circularized at radii ^ 1000 GM/c 2 . 

• The variability at energies ^ 5keV shows very low co- 
herence with the variability at energies ^ 5 keV. Given the 
the similarity of the PSD among the three energy bands, 
this low coherence is likely either due to a non-linear rela- 
tionship between the soft and hard energy bands, or due 
to multiple variability components some with strong corre- 
lations between soft and hard variability and other others 
with strong anti-correlations. 
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